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At low temperatures electron hopping in a three dimensional Coulomb glass produces fluctuations 
in the single particle density of states and hence in the resistivity. This results in a low frequency 
resisitivity noise spectrum which goes as f~ a where a is very close to 1. This holds down to 
extremely low frequencies. 
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Low frequency 1/f noise jl]-^ is ubiquitous; it is found 
in a wide variety of conducting systems such as metals, 
semiconductors, tunnel junctions and even supercon- 
ducting SQUIDs [|],||- Yet the microscopic mechanisms 
are still not well understood. In some cases the electrons 
are in a Coulomb glass which is an insulator with ran- 
domly placed electrons that have Coulomb interactions 
0. Lightly doped semiconductors and disordered met- 
als are examples of such systems. This paper focuses on 
1/f noise in Coulomb glasses. Experimental studies on 
doped silicon inversion layers have shown that low fre- 
quency 1/f noise is produced by hopping conduction || . 
More recent experiments have observed 1/f noise down 
to 0.1 Hz in boron-doped silicon ||. Because the sys- 
tems are glassy, electron hopping can occur on very long 
time scales and this produces low frequency noise. In this 
paper we show that the resulting noise spectrum goes as 
f~ a where / is frequency and the exponent a ~ 1. 

Shklovskii developed the first theory of 1/f noise in 
Coulomb glasses. He suggested that it is produced by 
fluctuations in the number of electrons in an infinite 
percolating cluster flQ|] . These fluctuations are caused 
by the slow exchange of electrons between the infinite 
conducting cluster and small isolated donor clusters. A 
more rigorous calculation combined with numerical simu- 
lations jll| of Shklovskii's model found a noise spectrum 
that went as f~ a where a was considerably lower than 1. 
Furthermore, below a minimum frequency of order 1-100 
Hz, the noise spectral density saturated and became a 
constant independent of frequency A similar conclusion 
holds for a model suggested by Kozub [ll| in which elec- 
tron hops within finite clusters produce fluctuations in 
the potential seen by hopping conduction electrons that 
contribute to the current. Kogan has argued that tran- 
sitions between valleys in the energy landscape produces 
1 /f noise because high barriers result in slow fluctuations 
in hopping conduction |13|. 

In our approach electron hopping shifts the single par- 
ticle energies e because they depend on Coulomb inter- 
actions with other sites. This leads to fluctuations in the 
single particle density of states g{e) which, in turn, pro- 
duces fluctuations in the conductivity. The conductivity 
depends on the density of states g(e w n) in the vicinity 



of the Fermi energy fi. Note that g(e « fi) can be af- 
fected by hops between sites i and j even if the energies 
on these sites are not near the Fermi energy because an 
electron or hole on site i or j can interact with other sites 
whose energy is (or was) near the Fermi energy. 

We start with a model of the Coulomb glass that fol- 
lows that of Baranovskii, Shklovskii, and Efros (BSE) 
]l4| . In this model, the electrons occupy the sites of a 
periodic lattice, and the number of electrons is half the 
number of sites. Each site has a random onsite energy <fn 
chosen from a uniform distribution extending from —A to 
A. Thus, <7 , the density of states without interactions, is 
flat. A site can contain or 1 electron. The Hamiltonian 
can be written as 



(i) 



where the occupation number m equals i if site i is oc- 
cupied and — ^ if site i is unoccupied, e is the electron 
charge, k is the dielectric constant and is the dis- 
tance between sites i and j. The single site energy is 
£i — 4>i + J2j ~£r~ n j- At zero temperature Coulomb inter- 
actions between localized electrons result in a so-called 
Coulomb gap in the single particle density of states that 
is centered at the Fermi energy |j] jij| , [l6] . 

We will use Mott's argument for variable range hop- 
ping [Q,|l7],|l8| to relate fluctuations in the density of 
states to fluctuations in the resistivity. One can regard a 
Coulomb glass as a random resistor network ]l9| with a 
transition between sites i and j associated with a resis- 
tance Rij given by 

Rij :=i?° exp(£y) (2) 

where the prefactor = fcT/(e 2 7? ) with 7^ given by 
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where D is the deformation potential, s is the speed of 
sound, d is the mass density, and A? — Ej — e% — e 2 / 'nrtj. 
is the change in energy that results from hopping from 
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i to j. ol — ko, b is the effective Bohr radius of a donor, 
and as is the usual Bohr radius (ag = Ti 2 /me 2 ). Wc 
will set the mass m equal to the electron mass so that 
as = 0.529A. In eq. (||), the exponent is given by 



Si: 



kT 



(4) 



The exponent reflects the thermally activated hopping 
rate between i and j as well as the wavefunction overlap 
between the sites. e,j is given by 0: 

13 y max[|ei - \ej - fj,\] , (si - n)(ej - fi) > 

(5) 

At both high and low compensations, electron hopping 
usually occurs on one side of the Fermi level p, and the 
lower expression applies. At intermediate compensations 
and in the regime of variable range hopping, hopping 
electrons often cross the Fermi level and the upper ex- 
pression applies. 

In the regime of variable range hopping Mott pointed 
out that hopping conduction at low temperatures comes 
from states near the Fermi energy. Let e = e — p,. If 
we consider states within e Q of the Fermi energy, then 
the concentration of states in this band is N(e ) = 
fl° 9{^)ds where g(e) is the density of states measured 
from the Fermi energy. So the typical separation be- 
tween sites is R = [iV(£ )] _1 / 3 . To estimate the re- 
sistance corresponding to hopping between two typical 
states of the band, we replace r%j with R and \sj — e,| 
with s in eqs. (||) and (||) to obtain £(e ). Minimiz- 
ing £(e ) yields e . Plugging this into eqs. (jj) and (||) 
yields the variable range hopping formula for the resis- 
tivity 7>(T) = p (T)exp(£(e )). 

In our model the noise results from electron hop- 
ping which produces fluctuations in the density of states 
g(e) = g(e) + Sg(e), where g(e) is the average density of 
states. This in turn creates fluctuations in N(e ), £(e )> 
e a , and p(T). We can calculate these fluctuations by ap- 
plying classical perturbation theory pFJ] i th e derivation 
of the variable range formula. To first order, <5£(e ) = 
5 P (T)/p(T) = -(2fcT ff (T,e ))- 1 fhjg(T,e)di. We 
have included the temperature dependence of the density 
of states because at finite temperatures the Coulomb gap 
fills in and the density of states no longer vanishes at the 
Fermi energy [pl|-p6[. The autocorrelation function for 
the fluctuations in the resistivity is 

<5p{TM)&p{T,t 1 )> 1 



p 2 (T) 4k 2 T*g 2 (T,e ) 
del de' < 5g{T,e,t 2 )5g(T,e',t 1 ) > (6) 
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We assume that there is no correlation between the fluc- 
tuations in the density of states at different energies, so 



< 5g(T, i, t 2 )Sg(T, t\) > = E < Sg(T, i, t2)Sg(T, e, t±) > 

6(e-e>) (7) 

where E is an energy of order 2e . Furthermore we as- 
sume that the time and energy dependence of the density 
of states autocorrelation function are separable, allowing 
us to write 



de < Sg(T, s, t 2 )Sg(T, e, h) >= C(e a , T)f(T, t 2 - h) 

(8) 



where we are assuming translational invariance in time 
(stationary processes). C(e 0l T) is a function of e and 
temperature. The function f(T,t) characterizes the time 
dependence of the return to equilibrium by the system 
after it is perturbed by a fluctuation in the density of 
states. Inserting eqns. (Q) and (||) in (||) yields 



<8p{T,t 2 )&p(T,t{) > 
p 2 (T) 



EC{e ,T) 
tfTV(T,£„) A ' 2 lJ 



(9) 



To relate this to the spectral density of the noise S(uj), 
let W(*2 - tx) =< 5p(T,t 2 )5p(T,ti) > and let ip p (co) be 
the Fourier transform of il>p{t 2 — According to the 
Wiener-Khintchine theorem || , for a stationary process 
the spectral density of fluctuations is given by p7[| 



S p (uj) = 2ip p (uj) 
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We do not know the temperature dependence of 
/(T, t), so for the moment we will suppress this and just 
refer to f(t). Theoretical calculations find that after 
large deviations from equilibrium, the density of states 
returns to equilibrium with a time dependence given 

by g(p,t) kit or g(ji,t) ~ t~ e where Q < 1 @. 

This agrees with experiments done at low temperatures 
p9| If we assume that these functional forms are 
also valid for f(t) which applies to small perturbations, 
then we obtain 1/f noise. We now describe the calcula- 
tion leading to this conclusion [p8[ . One starts with the 
Hamiltonian (|l|) but assumes that the Coulomb interac- 
tions arc turned on at time t = 0: 



H = } 4>im + > mn-jOU) 



(11) 



where the step function 9(t) is for t < and 1 for t > 0. 
So for t < the noninteracting density of states is a con- 
stant g a . Once the interactions are turned on, one follows 
the subsequent time development of the Coulomb gap. 
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The Coulomb gap arises because the stability of the 
ground state with respect to single electron hopping from 
an occupied site i to an unoccupied site j requires |t]] 
A\ > 0. So we need to subtract from the density of 
states those states which violate this stability condition. 
This leads to a self-consistent equation for the density of 
states @,||,|2[: 



g(e,t)=9oX\[l-al f 

j>i \ J - 



A 



di'g(i',t)6( + i - i') 



typical difference between successive iterations is typi- 
cally less than 1 part in 10 5 . We find that the Coulomb 
gap develops slowly over many decades in time p^| . Af- 
ter an infinite amount of time, the density of states at 
the Fermi energy /i goes to zero and g(i) ~ i 2 . 

The functional form of the time dependence of g(i 7 t) 
varies with the energy i and with g Q . For conduction 
noise we are interested in the time dependence of the den- 
sity of states at the Fermi energy which is shown in Figure 
[l] for 1CT 8 s < t < 10 8 s. For Jo = 2x 10 5 states/K-A 3 , 
we can fit our results to the form g(n,t) = £>i ln(i /i) 
where t < t , t a = 3 x 10 43 sec, and for g a = 6.25 x 10 5 
states/K-A 3 , g(fi, t) = B 2 t- e where 9 0.05. The val- 
ues of Bi, B2, and the other parameters used to obtain 
these results are given in the caption of Figure [I]. These 
fits change slightly for longer times. For example, for 
10" 8 s < t < 10 100 s and g a = 2 x 10 5 states/K-A 3 , 
the fit to our results has the form g(/J,,t) ~ t~ 6 where 
and n'j, respectively, while all other sites have their ^ 0.01. Still we see that the density of states at the 

Fermi energy approaches its equilibrium value roughly 
logarithmically in time. This is consistent with recent 
experiments on thin semiconducting 29 3^] and metal- 
lic films which have shown that the system adjusted 



Fin'i = 1, n' = O)0(t - m (e', i, r tJ )) 



where the single-site energy = i, ij = 



and 



(12) 



a n is 



the lattice constant, = rii + 1/2; so n ■ = 1 if site i is 
occupied and if site i is unoccupied. F(n£,r^) is the 
probability that donors i and j have occupation numbers 



ground state occupation numbers n' k . is the number 
of electrons which jump from site i to site j per unit time. 
6(t — Tij) represents the fact that at time t, the primary 
contributions to the change in the density of states will 
be from those hops for which Tij < t 33 1 . In writing eq. 
(|l2|), we assume that these hops together with phonons 
have equilibrated the system as much as is possible at 
time t. The hopping rate t^ 1 is given by |Q 

r^ 1 = 7 ° exp(-^)[l + N(Ai)]F(n' i = 1, n' 3 = 0) (13) 

where N(Aj) is the phonon occupation factor and reflects 
the contribution of phonon assisted hopping. We are also 
allowing for spontaneous emission of phonons since we are 
considering a nonequilibrium situation in which electrons 
hop in order to lower their energy. Following [^4 28 we 
can rewrite the self-consistent equation g(i,t): 



to changes in the Fermi energy approximately logarith- 
mically in time. These films were grown on insulating 
substrates which separated them from a gate electrode 
that regulated the electron density, and hence the chemi- 
cal potential, of the film. The conductance was measured 
as a function of the gate voltage. If the gate voltage was 
changed suddenly from, say, V to Vi, the conductance 
had a very fast initial rise, followed by a period of rapid 
relaxation, which in turn was followed by a long period of 
very slow relaxation. The relaxation could be described 
by Int or t~ e with 6 being small and varying slowly with 
time. This is consistent with our view |28| that when 



j-~ j A A di'g(e',t) £ drAnr 2 



*) = 9o exp < - 



F(n(s) = 1, n(e') = 0)6 (— + e - e')9(t - r(e', e, r)) } (14) 
nr J 

At low energies large distances play an important role 
and so we have replaced the sum by an integral over r in 
the exponent. The origin is at site i. n(e) is the occupa- 
tion probability of a site with energy e. r(e', s, r) is given 
by ( |l3| ) with replaced by r, £j replaced by e, and ij 
replaced by e' . 

Since it is not clear how the stability condition > 
can be applied to finite temperatures, we confine our cal- 
culations to the case of T = 0. In this case the phonon 
occupation factor N(Aj) = and the electron occupa- 
tion factor F(rii = 1, rij = 0) = 1 if ii < and ij > 0. 
Otherwise F(rii = l,rij = 0) = 0. We can solve eq. ( |l4| ) 
iteratively on the computer. After a few iterations the 



the gate voltage is changed, the Fermi energy changes, 
and time dependent relaxations arise because the system 
must dig a new Coulomb gap in the density of states at 
the new Fermi energy. 

So both theory and experiment indicate that the 
nonequilibrium density of states approaches its equilib- 
rium value roughly logarithmically in time. Returning to 
the original model described by (|l]) , we assume that this 
time dependence holds true in the linear response regime 
at low temperatures. If a fluctuation Sg(fi, t = 0) at t = 
pushes the density of states away from its mean equilib- 
rium value at the Fermi energy, then this perturbation 
will decay according to f(t) which enters into eqs. (||) 
and (^). Our nonequilibrium calculation indicates that 
f(t) can have the form: 



h{t)= B 1 ln(|) 



(15) 



where t < t a , and t a is on the order of the age of the 
universe or longer, or 
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h{t) = B 2 tr 



(16) 



where 9 « 1, and B\ and B 2 are positive constants of 
order g a . In both cases t is greater than some i m i n of 
order 10~ 8 s, say. The time dependence is a function of 
the energy, so here we set e = ji. Fourier transforming 
fi(t) and keeping the real part, we hnd that 



AM 



ttBx 1 



(17) 



This implies that the noise spectral density S(u>) ~ 1/lu. 
Fourier transforming f 2 (t) and keeping the real part 
yields 



AM 



ttB 2 8 1 



(18) 



for 6 < 1. This implies S'(w) - 1/w 1 - . 

To summarize, electron hopping leads to fluctuations 
in the density of states that relax back to equilibrium 
roughly logarithmically in time. This leads to 1/f noise 
in the spectral density S(u) of the noise in the resistivity. 
In particular we find that S(u>) ~ l/uj a where a = 1 if 
the relaxation is logarithmic in time, and a = 1 — 9 if the 
relaxation is a power law that goes as t~ e where 6 -C 1. 
In general a depends on temperature Q and is weakly 
dependent on the noninteracting density of states g and 
on the times scales. As eq. (|l0|) indicates, the noise am- 
plitude also depends on the temperature. Unfortunately 
we cannot ascertain these temperature dependences be- 
cause we do not know the temperature dependence of the 
fluctuations 6g(T,e,t) in the density of states. However 
we believe that our mechanism for 1/f noise should be 
valid at low temperatures (T ~ 20 K) where the loga- 
rithmic time dependence of the conductance is observed 
after the Coulomb glass has been pushed out of equilib- 
rium by the sudden application of a gate voltage [§9| [II) . 
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og =2x1CT 5 HC 1 A~ 3 
□ g =6.25x1(r 5 K~ 1 A~ 3 

1Cf 10 10- 5 10° 10 5 10'° 

time [s] 

FIG. 1. Density of states g(e = [i, t) at the Fermi energy as 
a function of time for different values of g - Solid lines are fits 
to the numerical integration of (^). The fit to the g a — 2x 10 5 
states/K-A 3 data is given by Bi ln(t /i) where lnt = 100 
with time measured in seconds and Bi — 1.64 x 10 _r 
states/K-A 3 . The fit to the g a = 6.25 x 10 5 states/K-A 3 
data is given by B 2 t~ e where 6 = 0.050 and B 2 = 1.82 x 10" 5 
states/K-A 3 . Parameters used are A = 2 x 10 4 K, k = 10, 
d = 7.18 g/cm 3 , s = 5.0 x 10 5 cm/sec, D = 5 x 10 3 K, and 
a = 4 A. The density d is chosen to be that of In203. 




